#include <iostream>
using namespace std;
const int N=2010,mod=1e9+7;
long long f[N][N];
int main()
{
    for(int i=1;i<=2000;i++)
    {
        f[i][1]=f[1][i]=i;
    }
    for(int i=2;i<=2000;i++)
    {
        for(int j=2;j<=2000;j++)
        {
            f[i][j]=(f[i-1][j]+f[i][j-1]+1)%mod;
        }
    }
    cout<<f[2000][2000]*f[2000][2000]%mod;

    return 0;
}

